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We investigate the scenario of homogeneous nucleation for a first order quark-hadron phase transi- 
tion in a rapidly expanding background of quark gluon plasma. Using an improved preexponential 
factor for homogeneous nucleation rate, we solve a set of coupled equations to study the hadroniza- 
tion and the hydrodynamical evolution of the matter. It is found that significant supercooling is 
possible before hadronization begins. This study also suggests that spinodal decomposition com- 
petes with nucleation and may provide an alternative mechanism for phase conversion particularly 
if the transition is strong enough and the medium is nonviscous. For weak enough transition, the 
phase conversion may still proceed via homogeneous nucleation. 
PACS number(s): 12.38.Mh, 64.60. Qb 



I. INTRODUCTION 

The hadronization of Quark Gluon Plasma (QGP) pos- 
sibly produced in the early universe or expected to be 
formed in relativistic heavy-ion collisions El has been 
the focus of much attention during the past few years. 
The quark gluon plasma (QGP), if formed, would ex- 
pand hydrodynamically and would cool down until it 
reaches a critical temperature Tc where a phase transi- 
tion from the quark matter to the hadron matter begins. 
Although the plasma has to hadronize, the mechanism 
of hadronization still remains an open question. The 
percolation model calculations are used in the case of 
a second order phase transition. In a first order scenario, 
the dynamics of the phase transition has been modeled 
in several ways. In the most idealized picture, the tem- 
perature of the plasma is held fixed at T = Tc until 
the phase conversion is completely over. Assuming an 
isoentropic expansion in (1+1) dimension, the hadroniza- 
tion in the above picture gets completed at t% = ttq 
where r is the ratio of the degrees of freedom of QGP 
and hadronic phases and tc is the proper time at which 
the QGP cools down to the temperature Tc- In real- 
ity, a first order phase transition is characterized by a 
large nucleation barrier that separates the two phases at 
T = Tc and the hadronization will not begin unless the 
matter supercools below Tc ■ Alternatively, the theory of 
homogeneous nucleation has been invoked to study the 
first order phase transition which is more realistic than 
the above idealized adiabatic scenario and has been in 
use for quite some time in the cosmological context [Q. 
In this picture, the transition is initiated by the nucle- 
ation of critical-size hadron bubbles from a supercooled 
metastable QGP phase. These hadron bubbles can grow 
against surface tension, converting the QGP phase into 
the hadron phase as the temperature drops below the 



critical temperature, Tc- For strong enough transition, 
the large amplitude fluctuations are suppressed so that 
the nucleation begins from a (nearly) homogeneous back- 
ground of supercooled metastable phase. This has been 
the basis of homogeneous nucleation theory j|] based on 
which the QCD phase transition has been studied ex- 
tensively [H— pi . However, for a weak enough transition, 
the matter may not remain in a pure homogeneous state 
even at T = Tc due to the transitions that may occur 
above Tc- Evidence of such pretransitional phenomena 
are quite common in several condensed matter studies 
particularly in the case of isotropic to nematic transition 
in liquid crystals. The possibility of such a transition 
had also been investigated during the cosmological elec- 
troweak phase transition in the early universe p0| , pl| . 
Due to high rate of thermal fluctuations and slow cool- 
ing of the universe, a strong phase mixing is expected 
even at temperature above Tc where the new phase is 
highly metastable. In such situation, the phase transi- 
tion either proceeds through percolation UM, if not, the 
dynamics of the first order transition is going to be quite 
different than the standard theory of homogeneous nu- 
cleation fl2|| . Similar phenomena is also expected in the 
case of a quark hadron phase transition both in the early 
universe as well as in the plasma expected to be formed 
during the relativistic heavy ion collisions [ p~3|JT^ ] . An 
ideal quark gluon plasma in (1+1) dimension expands as 
per the Bjorken scaling where T 3 r is constant p5| . This 
scaling would mean that the rate of change of tempera- 
ture is higher for the plasma at earlier time as compared 
to the plasma at later time. In the case of early uni- 
verse, due to high initial temperature, the rate of cooling 
becomes quite slow by the time it approaches Tc- How- 
ever, the QGP produced during the heavy ion collisions 
will have smaller initial temperature and will cool much 
more rapidly as compared to the early universe. An- 
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other difference as compared to the early universe is that 
the QGP produced at RHIC and LHC energies may at- 
tain kinetic equilibrium in a very short time « 0.3 — 0.7 
fm/c but will remain far off from chemical equilibrium 
Jjq , [j"7| . Such a plasma will be more gluon rich and many 
more quark and anti-quark pairs will be needed before 
the plasma achieves chemical equilibrium. The chemi- 
cally unsaturated plasma will cool still at a faster rate 
since additional energy will be consumed in approaching 
the chemical equilibration (ljjl^]. Recently, we have in- 
vestigated the effect of thermal fluctuations leading to 
phase mixing by modeling them as subcritical hadron 
bubbles p9[ . Although the equilibrium density distri- 
bution of these subcritical bubbles can be quite large, 
their equilibration time scale is larger than the cooling 
time scale for the QGP. As a consequence, for RHIC and 
LHC energies, they will not build up to a level capable of 
modifying the predictions from homogeneous nucleation 
theory: the QGP has to supercool and the hadronization 
may still proceed through the nucleation of critical bub- 
bles within a nearly homogeneous background of quark- 
gluon plasma. 

For the QGP formed at RHIC and LHC energies, the 
question to be asked is how long the plasma would con- 
tinue to supercool. It is expected that the plasma will 
supercool upto some temperature T m until the density of 
nucleated hadron bubbles become sufficient to heat up 
the medium due to the release of latent heat Since 
the medium is heated up, the nucleation is again switched 
off at some point and the hadronization continues due to 
the growth of previously nucleated hadron bubbles. Re- 
lated to this another crucial aspect which needs investiga- 
tion is that how fast the barrier between the metastable 
QGP phase and the stable hadron phase decreases as 
the system supercools. For strong enough supercooling, 
the barrier between the stable and metastable minima 
completely vanishes leading to a point of inflection at 
T = Ts known as spinodal instability. Thus, the rapidly 
quenched system leaves the region of metastability and 
enters the highly unstable spinodal region before a sub- 
stantial amount of nucleation begins. The spinodal de- 
composition has also been suggested as a possible mecha- 
nism of phase conversion for a rapidly expanding system 
of quarks and gluons p(J . The second scenario of phase 
transition may lead to coherent pion emission due to the 
formation of disoriented chiral condensates pl[ . 

For the case of QGP produced in heavy ion collision, 
the route to hadronization either through nucleation or 
through spinodal decomposition depends sensitively on 
several factors like nucleation rate, effective potential 
used to model the dynamics of the phase transition and 
also on the hydrodynamical evolution that decides the ex- 
pansion rate of the matter. The homogeneous nucleation 
rate is estimated from the factor I — Iq exp(— AFc/T) 
where Iq is the product of dynamical and statistical pref- 
actors and AFq is the minimum energy needed to cre- 
ate a critical hadron bubble. The exponential factor 
dominates the rate at temperature close to Tq whereas 



the prefactor Iq plays the significant role away from Tq. 
The most commonly used expression for the prefactor is 
Io ~ T 4 lEi which may be alright for the case of early 
universe but not for the heavy ion collisions where we 
expect large amount of supercooling. The theory of ho- 
mogeneous nucleation assumed significance for the study 
of phase transition in QGP produced in the Rclativistic 
heavy ion collision after Csernai and Kapusta j| derived 
an improved prefactor using coarse grained field theory. 
An important aspect of their formalism is the presence of 
viscosity which is essential for the removal of latent heat 
away from the hadron bubble. The prefactor vanishes 
for an ideal plasma. This would suggest that a dissipa- 
tive dynamics should be followed for the phase transition 
to be completed through nucleation In a subsequent 
work j^], Ruggeri and Friedman, on the other hand, de- 
rived a prefactor which does not vanish in the limit of 
zero viscosity. Recently, we have also derived a prefactor 
which is more general and also reproduces the Csernai- 
Kapusta and Ruggcri-Fricdman results in the limiting 
cases |23| . The prefactor Iq in all these estimates is found 
to be less as compared to T A . In this work, we estimate 
the nucleation rate using an improved prefactor as de- 
rived in Ref. pi. The argument AFc/T that appears 
in the exponential depends on the phenomenological po- 
tential model used to describe a first order phase tran- 
sition. Since the lattice predictions |M 25 are not yet 
conclusive about the nature and the strength of the tran- 
sition, we use an effective potential that follows the bag 
equation of state and also covers a wide range from very 
weak to very strong first order phase transitions. Finally, 
we study amount of supercooling and rate of hadroniza- 
tion by solving self consistently the nucleation rate along 
with energy momentum conserving hydrodynamic equa- 
tion both for dissipative and nondissipative plasma. The 
dissipative plasma makes the evolution slow and also gen- 
erates extra entropy. From this study we show that spin- 
odal decomposition may compete with homogeneous nu- 
cleation if the plasma is nondissipative and the transition 
is relatively stronger. Since the amount of supercooling 
is less in case of a dissipative plasma, the phase conver- 
sion may still proceed through the nucleation of critical 
size hadron bubbles from a homogeneous background of 
supercooled quark gluon plasma. 

The paper is organized as follows: In Sec. II, we briefly 
review the bag equation of state and an effective poten- 
tial used to model the phase transition. The nucleation 
rate with various prefactors are discussed in Sec. III. Fi- 
nally, we present our results in Sec. IV followed by the 
conclusion in Sec. V. 



II. THE EFFECTIVE POTENTIAL AND THE 
EQUATION OF STATE 

The spinodal instability corresponds to the inflection 
point in the effective potential (below Tc) used to study 
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the dynamics of a phase transition. Obviously, the tem- 
perature Ts at which the instability occurs will depend 
on the parameters of the effective potential. Since the or- 
der as well as strength of the quark hadron phase transi- 
tion is still an unsettled issue, we consider a more generic 
form of the potential which covers a wide range from very 
strong to very weak first order phase transition. The pa- 
rameters of the effective potential are determined from 
the physical observables like the surface tension a, the 
correlation length £ and also from the requirement that 
the potential difference between the two minima should 
corresponds to the pressure difference between the quark 
and hadron phases. We follow the bag equation of state 
to estimate various thermodynamical observables. For 
example, the pressure can be estimated from the relation 
P = TY,\nZi/V — B where the partition function In Zi 
for a single fermion or boson is given by |^7j] 
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ln(l±e-V^?/ T )] , (1) 



where the + and — stand for fermions and bosons respec- 
tively. In the above, the chemical potential \x has been 
set to zero. Assuming, a baryon free quark gluon plasma 
which consists of it, d, and s quarks and gluons, the total 
pressure can be written as 
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(2) 



where n/ is the number of flavours and rtii, (i=u, 
d, s) are their masses. The function /(m^/T) = 

(360/7tt 4 ) J °° du it 2 In (l + e -V^+KWj and B is the 

bag constant. In Eq. (^), the first term is due to massless 
gluons while the second term needs to be evaluated for 
different physical masses of the quarks. If the quarks are 
taken as massless, we get 
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If the physical masses of quarks are used as m u ~ ~ 8 
MeV and m s ~ 160 MeV, then at the critical tem- 
perature (T c = 160 MeV) we get f{m u /T) ~ 1 and 
f(m s /T) ~ 0.8. Thus, 
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which is not much different from the one for the massless 
gas. Similarly, the total pressure in hadron phase con- 
sisting of n different types of particles and resonances is 
given by 
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Here. 

h(mi/T) 



.« J- d MU 2 m(l- e -\/" 2 +(™.AT). If we 

consider only a massless pion gas, the above summation 
contributes only a factor of 3. Assuming the hadron gas 
consisting of 7T, K, tj, p and u) with appropriate mass 
and degeneracy [m^ = 137.0 MeV (g=3), rriK — 496.0 
MeV (g=4), m n = 547.45 MeV (g=l), m p = 768.50 MeV 
(g=9), m w = 781.94 MeV (g=3)] we get 

TT 2 

p h ~ a h T 4 , where a h = 4.5— at T = I60MeV. (6) 

90 

In addition to the above equation of state, we use 
the following phenomenological potential to describe the 
phase transition p6j 



V(<fi, T) = a(T) (f) 2 - bT(f> 3 +c(t) 4 , 



(J) 



where b and c are positive constants. Although, the ex- 
act knowledge of the order parameter <f> is not necessary 
for the present purpose, it can be related to entropy or 
energy density Q or to the sigma and pion field |2(J . The 
above potential has two minima, one at <j) q = and the 
other one at 



(36T + v / 9& 2 T 2 -32ac)/8c 



(8) 



which in the present case will represent quark and hadron 
phases respectively. These phases are separated by a 
maximum that occurs at <f) m given by 



m = (36T - y/9b 2 T 2 - 32ac)/8c. 



(9) 



In the thin wall approximation |22| , b and c can be ex- 
pressed in terms of surface tension a and correlation 
length i as pj 



b = 



1 



12f 3 CT 



(10) 



The height of the degenerate barrier at T = Tc or at 
a(T c ) = b 2 T c /Ac is given by 



V b (T c 
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(11) 



Thus, the potential can describe first-order phase transi- 
tions, with varying strength by either changing a or £ or 
both. The latent heat is also a measure of the strength 
of the first order phase transition which in the bag model 
is given by AB. The requirement that at all the temper- 
atures, the difference between the two minima should be 
equal to the pressure difference between the two phases 
Q fixes the third parameter a(T) as 
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Ap= Ph - Pq = -[v(<j> h )-v(a)] 

= -[a(T)-W> fc + c$] eg. 



(12) 



Scavenius and Dumitru pO| have used a linear sigma 
model (LSM) potential and fitted the potential of the 
form given by Eq. ([?]) to fix the parameters a, b and c. 
With this form of potential they obtain surface tension a 
and correlation length £. We directly use surface tension 
(a) and correlation length (£) to fix the parameters b and 
c. For any value of b and c the nonperturbative vacuum 
effect is taken care of by bag constant B. At T = Tc, 
$ h = bT c /2c = y/Gag. For a = 10 MeV/fm 2 and £ = 0.7 
fm (typical hadron size) we get 0/, = 91 MeV which is 
nothing but the pion decay constant. 

Figure 1 shows the plot of V(4>) as a function of 4> a t 
four different temperatures for typical values of a = 30 
MeV/fm 2 and £(T C ) = 0.7 fm. At T = T c , the poten- 
tial is degenerate with V(4> q ) = V(<ph) = being sepa- 
rated by a barrier at (f> = <p m . As the matter supercools, 
the hadron phase has lower free energy as compared to 
the QGP phase which is held fixed at <f> q = and also 
V(4> q ) = 0. Thus, below Tc, the cf> q phase is metastablc 
with respect to the stable hadron phase at <f> = cj>h. At 
T = Ts, the potential develops an inflection point where 
4> q and 4>m coincide. The condition <p q — cf> m = leads to 
a(Ts) = and <f>h — 36T,g/4c . The spinodal tem per ature 
Ts can be obtained analytically by solving Eq. ( |l2| ) as 
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At T = Ts their exists only one minimum correspond- 
ing to the hadron phase. If the QGP supercools upto 
this point it will become unstable and may go to hadron 
phase by spinodal decomposition. It is worth noting that 
as tr — ► and £ — > oo, Ts — > Tc- The spinodal temper- 
ature depends on the strength of the transition. For a 
strong enough transition, a is large and Ts is lower as 
compared to the case when the transition is weak. We 
are interested to know whether the system cools down 
to the temperature T5. For comparison we denote the 
minimum temperature reached during the supercooling 
as T m . It is the temperature in the supercooling region 
at which the system starts reheating due to the release of 
latent heat. The rate of nucleation will be suppressed for 
a stronger first order phase transition due to large bar- 
rier resulting in higher supercooling, (i.e., smaller value 
of T m ). Thus, both T5 and T m depend on the strength of 
the transition and need to be evaluated properly. While 
Ts can be estimated directly from Eq. ( |l3| ) , T m requires a 
self consistent solution of a set of equations involving the 
nucleation rate and energy momentum conserving hydro- 
dynamical equations. 
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FIG. 1. The free energy density V{4>) as a function of <j> at 
four different temperatures. 



III. NUCLEATION RATE 

The homogeneous nucleation theory assumes the for- 
mation of nucleating clusters within the initially homo- 
geneous metastable state. One can use the Langer for- 
malism ||, to calculate the probability of nucleation of 
hadron bubbles per unit time per unit volume from a 
homogeneous background of quark gluon plasma, given 
by 



I = I Q exp 



AF C 
T 



(14) 



Here, AFc is the minimum energy needed to form a crit- 
ical hadron bubble and Iq is the prefactor which can be 
written as Iq = Kf2/27r. The statistical factor is a mea- 
sure of both the available phase space as the system goes 
over the saddle and of the statistical fluctuations at the 
saddle relative to the equilibrium states. The dynami- 
cal prefactor k gives the exponential growth rate of the 
bubble or droplet sitting on the saddle. 



A. Statistical prefactor 

To understand the meaning of the statistical pref- 
actor, consider a classical system with N degrees of 
freedom described by a set of N collective coordinates 
rji, i — 1, . . . , N. The coarse-grained free energy func- 
tional F{r/} of the system has local minima F{rji} in the 
{?/}-space, corresponding to the metastable {77 } and sta- 
ble states , separated by the energy barrier. The point 
of minimal energy along the barrier is the so-called sad- 
dle point {rj }. The rate of the decay of the metastable 
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state is determined by the steady-state current across the 
saddle point from the metastable to the stable minimum 
of F {r]}. According to ]3|, p8| , p9| , the statistical prefactor 
can be written as 



»7T(T 

3IAJ 



3/2 



/2ttT 



1/2 



det(M /27rT) 



det{Mi/2TiT) 



(15) 



where Ai is the negative eigen value of the matrix My = 
d 2 F / drj i dr]j 1 evaluated at the saddle point. The sub- 
script denotes the metastable state and the prime in- 
dicates that the negative eigen value Ai as well as the 
zero eigen values of the matrix My is omitted. The cal- 
culation of the fluctuation determinant in Eq. (|l^) is ex- 
tremely difficult and adds to the uncertainty of the cal- 
culation of SI. In Ref. pHl, the above prefactor has been 
estimated under harmonic approximation using the prod- 
uct of eigen values of the mobility matrix M, evaluated 
at the saddle point and the metastable points. The final 
expression that accounts for fluctuation correction reads 



3\/3 



3/2 



Rc 



(16) 



where a is the surface tension, £ is the correlation length 
and Rc is the size of the critical droplet. Since Rc is in- 
finite at T = Tc, Sli diverges as T — > Tc- Alternatively, 
one can follow the earlier approach of Langer || where 
the fluctuation corrections are absorbed into the free en- 
ergies of the metastable (F ) and saddle point region (Fs) 
and the activation energy of the critical droplet is sim- 
ply AFc — Fs — Fo. Therefore, omitting the fluctuation 
determinant, the expression for is given by 



327T 2 r 1 / 2 /CT\3/2 



(17) 



In the above, a is interpreted everywhere as the true sur- 
face energy that includes fluctuation corrections. This 
approach has been adopted in Refs. to estimate the 
statistical prefactor. In order to evaluate Ai, we will ap- 
proximate the potential barrier between the metastable 
and stable states by the excess of the free energy AF 
corresponding to the formation of a spherical bubble of 
radius R. In the thin wall approximation |2lq ], it can be 
written as the sum of the bulk and the surface energies 



AF(R,T) 



47T 



-R 3 Ap + 4irR 2 <r, 



(18) 



where Ap is the pressure difference between the hadron 
and QGP phases. Minimization of AF with respect to 
the radius R yields the free energy of the critical bubble 



AFr 



-7ri? c er; R c = 



2a_ 
Ap' 



(19) 



It is convenient to introduce a similarity number A 2 
Rcx/iira/T and the reduced radius r = R/R c @,|| 



Two functions AFi(i?) and Ai^(-R) are similar if they 
have same similarity number. In the harmonic approx- 
imation, Eq. (U|) can be expanded around the critical 
radius as 



AF _ (AF\ | 1 / d 2 AF 

~ \~)r=r c + ^ R=Rc 

= ^-A2(r-l)' 



{R-Rcf 
(20) 



and we get finally for the negative eigen value Ai 

Ai = -2TA 2 . (21) 

Following [0J|], we estimate using Eq. @ and 
Eq. ©. 

B. Dynamical prefactor 

Csernai and Kapusta Q have derived a dynamical pref- 
actor k by generalizing the coarse-grained effective field 
theory of Langer to a relativistic case. For a baryon-free 
plasma, where the thermal conductivity vanishes, the dy- 
namical prefactor is found to depend on the viscosity co- 
efficients and is given by 



kck 



4(7(4^/3 + 
(Acj) 2 R c 



(22) 



Here r\ q and £ g are the shear and bulk viscosity coeffi- 
cients of the quark phase respectively and the Au> is the 
difference between the enthalpy densities of the QGP and 
the hadron phases, u> = e +p. According to Csernai and 
Kapusta (CK), there will be no bubble growth in the 
case of an ideal plasma with zero viscosity. Since vis- 
cosity is an essential ingredient in the above expression, 
for consistency, dissipative hydrodynamics should also be 
used to describe the space time evolution of the matter 
H . The CK approach has also been extended to include 
thermal dissipation in addition to viscous damping for 
the case of baryon rich quark gluon plasma |j2j . On the 
contrary, Ruggeri and Friedman || argued that for rela- 
tivistic hydrodynamics the energy flow does not vanish in 
the absence of any heat conduction or viscous damping. 
Since the change of the energy density e in time is given 
in the low velocity limit by the conservation equation 
0, de/dt — -V.(wv) which implies that the energy flow 
oc low is always present. Therefore, following a different 
approach, Ruggeri and Friedman derived an expression 
for k, given by (for zero viscosity) 



krf = 



2<7LU„ 



R%(Alj) 



1/2 



(23) 



The above result is in contradiction with the expression 
given by Eq. (22). Recently, following the Langer's pro- 
cedure, we have solved the relativistic hydrodynamic in 
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the hadron, quark and the interfacial regions to obtain 
the dynamical prefactor. The k as obtained in Ref. [f23| , 
is given by 



k = Kq + k v : Kq = c s 



3R C 



^ 1 f*k +c 



6i? c UJ q 



(24) 



Here c s is the velocity of sound in the massless gas. Under 
certain assumption for c s , the above prefactor can also 
be written as |23| 



2auj„ 



4N ! 



1/2 



1 a(4V3 + C g 
c 2 (A^R C 



(25) 



As can be seen, the first term in the above equation is the 
same as obtained by Ruggeri and Friedman correspond- 
ing to the case of nonviscous plasma [see Eq. (p3|)]. The 
second term is similar to the result obtained by Csernai 
and Kapusta. Two important aspects of our result are 
(a) the prefactor k can be written as a linear sum of a 
nonviscous (kq) and a viscous (k„) components and (b) 
the nonviscous component (kq) which depends on two 
parameters R and £ is finite in the limit of zero viscos- 
ity. Further, it has been argued in Ref. [^3| that the 
nonviscous part of the prefactor basically arises due to 
nonuniform pressure across the interface. In the present 
work, we will use Eq. ( p4| ) to calculate the dynamical 
prefactor. 



IV. NUCLEATION AND THE EXPANSION 
DYNAMICS 

To calculate the fraction of space converted to hadron 
phase using the nucleation rate I(T), one requires a ki- 
netic equation. The simplest of such equation is given in 
Ref. J6J. According to it, if the QGP cools to Tc at a 
proper time tc, then at some later time r the fraction 
h of space which has been converted to hadronic gas is 
given by 



h(r) 



dT'I(T(r'))[l-h(r')}V(r',T). 



(26) 



Here V(r', r) is the volume of a bubble at time r which 
had been nucleated at an earlier time r'; this takes into 
account the bubble growth. The factor [1 — /i(r')] is the 
available space for new bubbles to nucleate. The model 
for bubble growth can be taken as ||J3^] 

V(r', r) = y U(T(t')) + £ dr"v{T{r"))\ " , (27) 

where v(T) is the velocity of the bubble growth at tem- 
perature T and is taken to be v(T) = 3[1 - T/T c } 3/2 
P,p3. This expression is intended to apply only when 



T > \Tc so that the growth velocity stays below the 
speed of sound of a massless gas, c/ At the critical 
temperature, both the nucleation rate and growth rate 
vanish. One can also write v in terms of pressure differ- 
ence between the two phases as used in Ref. [[}5| . 

The evolution of the energy density is given by 



15p 37 1 



de Dlo _ Ar/f'S + C _ A* 
dr t t 2 t 2 



(28) 



where D—l for the expansion in (1 + 1) dimension. The 
factors r) and ( are the shear and the bulk viscosity of 
the medium. For nonviscous plasma (for zero viscosity), 
the above equation follows the Bjorken's scaling solu- 
tion where T 3 r is a constant. The energy momentum 
equation needs to be solved numerically for expansion in 
(3 + 1) dimensions. However, retaining the simplicity, we 
can still use Eq. ( p8| ) for spherical expansion |3|^] with 
the choice of D = 3, although the scaling u M = /r 
(where u M is the four velocity of the fluid) may not be 
valid if viscosity is included. 

In the transition region, the energy density at a time 
t can be written in terms of hadronic fraction h(r) as 



;(r) = e,(T) + [e h (T) - e q (T)]h(r). 



(29) 



Here e q = 3p q + AB and eu = Sph are energy densities in 
the QGP and hadron phase respectively. The viscosity 
at any time r can also be written pZfl as 



/i(r) = fi g (T) + \fi h (X) - » g (T)]h(T), 



(30) 



where \i q = ii] q /3 + ( q and fj,h — ^Vh/i + Ch- The tem- 
perature dependence of the viscosities are taken as ~ T 3 

M. 



V. RESULTS AND DISCUSSIONS 

We compare the nucleation rates using different dy- 
namical prefactors k as given by 



il 2 KcK ( AF C 
i\ = — exp 



h = 

h = 



2?r 
2tt 
2tt 



exp 



exp 



T 
AF C 

T 
AF C 
T 



fJ 2 (K + Kv) 



2ir 

h = T 4 exp ( - 



exp 



AFc 
T 



AF C 
T 



(31) 



In all the cases, we use the same statistical prefactor 
which is more physical except for I5 where the to- 
tal prefactor is of the order of ~ T 4 as commonly used 
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in the literature ||22(] . Further, we choose £ = 0.7 fm and 
use Eq. ( |l9| ) to estimate AFc/T. The nucleation rates /2 
and Is correspond to the case of a nondissipative plasma 
whereas viscosity enters as an essential ingredient for the 
evaluation of I\ , I3 and I4. Figures 2 and 3 show the plot 
of above nucleation rates as a function oiT/Tc for a = 10 
MeV/fm 2 and a — 30 MeV/fm 2 respectively each at two 
typical values of viscosity coefficient r\ q . We consider only 
the shear viscosity r\ q since the bulk viscosity ( q that pro- 
vides resistance to expansion (or contraction) does not 
exist in the case of an incompressible fluid. Further, the 
viscosity coefficients rjh and C/i m hadron phase are ne- 
glected since the amount of supercooling is affected by 
viscosity of quark phase only. Notice that the most com- 
monly used prefactor ~ T 4 overpredicts the nucleation 
rate (J5) over a wide range of temperature particularly 
for small amount of supercooling and also when the tran- 
sition is relatively stronger (see Fig. 3). Further notice, 
that both for strong and weak transitions but for mod- 
erately viscous medium, the nucleation rates are domi- 
nated by the nonviscous components (since I4 is nearly 
the same as I 2 , see Figs 2(a) and 3(a) for rj q = 5T 3 ). The 
rates 1\ and I3 which contain viscosity have only a sec- 
ond order effect unless the medium is highly dissipative. 
At ?] q = 14. 4T 3 , the viscous component J3 competes with 
the nonviscous component I2 although 1\ is still lower, 
[see Figs. 2(b) and 3(b)]. From the above studies we may 
conclude that the choice of T 4 overestimates whereas the 
rate with the Csernai-Kapusta prefactor underesti- 
mates the nucleation rates in the region of interest as 
compared to 1 4. Since the nucleation rate J4 has both 
dissipative and nondissipative components, we will use it 
to estimate the nucleation rate and supercooling in the 
subsequent study. 
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FIG. 2. Nucleation rates obtained using different prefactors 
for two viscosities using surface tension a = 10 MeV/fm 2 . 
The expressions for Ii etc. are given in the text. 
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FIG. 3. Same as Fig. 2 but with a = 30 MeV/fm 2 . 

Next, we solve the coupled equations Eq. ( |26| ) and 
Eq. (g8[) to study the space time evolution of the mat- 
ter. The critical temperature is fixed at Tc = 160 MeV 
which gives i? 1 / 4 = 236 MeV. The strength of the tran- 
sition depends on a and £ or more explicitly on the ratio 
a 7£. Recent lattice calculations |2j| predict a between 
2 to 10 MeV/fm 2 . These calculations are without any 
dynamical quarks and thus are only indicative. There- 
fore a is treated as a free parameter in the present study 
and £ = 0.7 fm seems the most reasonable value for the 
correlation length. 

Another important parameter that affects the evolu- 
tion is the time tc that the plasma takes to cool down 
to T = Tc- Obviously, tq will depend on the initial 
temperature (X^), the formation time (r^) as well as on 
the expansion dynamics of the plasma. In the Bjorken 
scenario tq — (Ti/Tc) 3 Ti. It is quite reasonable to 
assume that the plasma will expand in (1 + 1) dimen- 
sion between r, and tc until the longitudinal dimension 
becomes comparable to the size of the colliding nuclei. 
Assuming Tj = 320 MeV and Tj = 1 fm/c, the Bjorken 
scaling predicts tc = 8 fm/c for expansion in (1 + 1) di- 
mension. As mentioned before, tq is a crucial parameter 
that affect the solution of the coupled equations which 
in turn depends on the initial conditions as well as on 
the expansion scenario of the plasma. Since tq is not 
known exactly, we assume it of the order of 6 to 8 fm/c. 
From tq onwards, we consider the spherical expansion. 
This expansion scenario corresponds to the fastest cool- 
ing rate resulting in maximum supercooling. Therefore 
the choice of D = 3 provides a lower bound on T m which 
we compare with the spinodal temperature Tg. 

Figure 4 shows the plot of T/Tc as a function of r for 
a nonviscous plasma at a typical value of tc = 8 fm/c. 
Unlike the ideal Maxwell construction ( where the tem- 
perature is held fixed at T = Tc until phase conversion 
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is over), the system supercools upto T m . At T m , the 
number density of the nucleated hadron bubbles is suf- 
ficient to raise the temperature again due to the release 
of the latent heat. The amount of supercooling will ob- 
viously depend on the strength of the transition being 
more for a stronger transition (large value of a) as com- 
pared to that for the weaker one (small value of a). A 
faster expansion also reduces the number density of the 
nucleated hadron bubbles resulting in a larger supercool- 
ing. Although shown for tc — 8 im/c, the amount of 
supercooling will increase further if tq is reduced due to 
faster expansion. We have also done the calculations for 
tc = 6 fm/c which are summarized in Fig. 6. For tc still 
smaller the expansion in (1+1) should be used. 
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FIG. 4. The temperature evolution obtained with the nu- 
cleation rate I4 for a spherically expanding system for differ- 
ent surface tensions. 

Next, we consider the dissipative hydrodynamics. As 
mentioned before, the viscous contribution to the nucle- 
ation rate is not very significant if the transition is strong 
enough and also the medium is moderately viscous. How- 
ever, the presence of a small amount of viscosity will 
affect the hydrodynamic evolution of the plasma. The 
temperature of the plasma will fall at a slower rate due 
to viscous heating of the medium. Figure 5 shows a plot 
of T/T c as a function of r for ^=0, 5T 3 and 14.4T 3 . 
Here we have considered the case of a relatively stronger 
transition (a = 30 MeV/fm 2 ) where the viscosity con- 
tribution to the nucleation rate is not very significant. 
However, significant effect can be seen on T m which in- 
creases (less supercooling) with increasing viscosity as 
expected ||. 
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FIG. 5. The temperature evolution obtained with the nu- 
cleation rate I4 for a spherically expanding system for differ- 
ent viscosity coeffs. 

Figures 6 and 7 show the plot of T m /Tc as a function 
of o at Tc= 6 fm/c and 8 fm/c respectively both for 
ideal as well as for dissipative plasma. The amount of 
supercooling becomes less as the medium becomes more 
viscous. The spinodal temperature T5 which depends on 
the ratio <r/£ has also been plotted in the same figures 
for comparison. First consider the case of nondissipative 
plasma with zero viscosity. The curves T m and T$ show 
a cross over point depending on the choice of tc- The 
nucleation rate is suppressed for a strong enough transi- 
tion due to large nucleation barrier resulting in a higher 
amount of supercooling. On the other hand for weak 
enough transition, the amount of supercooling is smaller 
and also well above the point of spinodal decomposition. 
Since T m depends on expansion rate of the medium, the 
cross over point will sensitively depend on tc] moving 
towards left for the faster expansion. Qualitatively it 
can be concluded here that the homogeneous nucleation 
is still the dominant mechanism of phase conversion if 
the transition is weak enough. For stronger transition 
and fast enough cooling, the system may reach T5 and 
the phase conversion may proceed through the spinodal 
decomposition. If the plasma is viscous, the amount of 
supercooling is reduced further due to slow evolution of 
the medium. Even the cross over point also shifts to- 
wards right showing that the nucleation is the dominant 
mechanism over a wide range of a /£ ratios. Due to un- 
certainties in tc, <r/£ and r\ q , it is difficult to say what is 
the exact scenario at RHIC and LHC energies. Assuming 
a weakly first order transition for a in the range of 2 to 
5 MeV/fm 2 , £ - 0.7 fm ||, 77, - 2T 3 (lower limit glj) 
and tc ~ 4 to 8 fm/c |lS| ], the nucleation still seems to be 
the dominant mechanism of phase conversion as opposed 
to the spinodal decomposition. 
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It may be mentioned here that we have considered a 
dynamical growth rate which is exponential in nature. 
This may not be true at the later stage of the expansion. 
The fusion of bubbles may also play a role at the later 
stage when the density of bubbles is sufficiently high [B8| . 
However, in the study of supercooling, we are very much 
in the domain of exponential region where the effect of 
bubble fusion can be ignored. 
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VI. CONCLUSIONS 

We have investigated the mechanism of phase conver- 
sion from quark gluon plasma phase to hadron phase via 
two routes: the standard homogeneous nucleation and 
spinodal decomposition. The point of spinodal decompo- 
sition depends on the strength of the transition or more 
precisely on the ratio a /£ where a and £ are the surface 
tension and correlation length respectively at T = Tc- 
The nucleation and supercooling, on the other hand, de- 
pend on the strength of the transition as well as on the 
expansion dynamics of the medium. We have solved a 
set of coupled equations to estimate the amount of max- 
imum supercooling under spherical expansion scenario. 
Which way the hadronization will proceed depends sen- 
sitively on the nucleation and expansion dynamics since 
the nucleation and expansion time scales are comparable. 
Qualitatively we can describe the results as follows: For 
strong enough transition with zero or very small amount 
of viscosity, the system reaches the spinodal instability 
before the amount of nucleated hadron bubbles become 
significant to begin phase conversion. The phase con- 
version in such a case will proceed via spinodal decom- 
position. If the medium is viscous or the transition is 
weak enough or both, the supercooling is much less; The 
phase conversion may still proceed through homogeneous 
nucleation. However, depending on the range of the pa- 
rameters, there could be a competition between the ho- 
mogeneous nucleation and the spinodal decomposition. 
A definite answer can be provided only when the param- 
eters such as surface tension, tc, viscosity and also the 
expansion scenario are known precisely. 
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